Entropy of solid 4 He: the possible role of a dislocation glass 
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Solid 4 He is viewed as a nearly perfect Debye solid. Yet, recent calorimetry indicates that its 
low-temperature specific heat has both cubic and linear contributions. These features appear in 
the same temperature range (T ~ 200 niK) where measurements of the torsional oscillator period 
suggest a supersolid transition. We analyze the specific heat to compare the measured with the 
estimated entropy for a proposed supersolid transition with 1% superfluid fraction. We find that 
the experimental entropy is substantially less than the calculated entropy. We suggest that the low- 
temperature linear term in the specific heat is due to a glassy state that develops at low temperatures 
and is caused by a distribution of tunneling systems in the crystal. It is proposed that small scale 
dislocation loops produce those tunneling systems. We argue that the reported mass decoupling is 
consistent with an increase in the oscillator frequency as expected for a glass-like transition. 

PACS numbers: 73.21.-b 



Introduction. A supersolid is a novel state of matter 
that simultaneously displays both superfluidity and crys- 
talline order. 4 Hc is thought to be a most likely candi- 
date for the supersolid state. Recent torsional oscillator 
experiments by Kim and Chan Q, Q generated renewed 
interest in this possibility. Pioneering work by Andreev 
and Lifshitz 0, Reatto p|, Chester [jj, Leggett 0, and 
Anderson Q laid the foundation for our thinking about 
this enigmatic state. Recently, Anderson et al. revisited 
this problem 0,0] and latest developments, presented at 
a KITP workshop, are available online 01 ■ 

In addition to the work of the PSU group 0,0, there 
is now an independent confirmation of the anomalous be- 
havior of solid 4 He, as presented by the groups of Reppy 
and Shir ahama [ill fl2T] . Both groups use torsional oscil- 
lators similar to the one by the PSU group of Chan 0,0] 
At the same time Rittner and Reppy [11| reported a his- 
tory dependence of the signal, when annealing the sam- 
ple, to the extent of no observation of any mass decou- 
pling in the torsional oscillator experiment. These tor- 
sional oscillator experiments clearly indicate anomalous 
mechanical properties of solid 4 He at low temperatures. 
However, the relationship between the mechanical mea- 
surements and the suggested superfluidity is not direct. 
The most direct proof of superfluidity would be obser- 
vation of persistent current. In this regard, we mention 
a recent experimental search for superflow by Beamish 
et. al. ^| that indicates no mass flow of any kind to 
very high accuracy. Thus, we conclude that the effect 
first observed by the PSU group is likely not an intrinsic 
property of solid 4 He, because it depends critically on the 
3 He concentration and shows annealing dependence pdl |. 

Alternatively, we suggest that many of the experimen- 
tal facts seem to be consistent with a glass-like behavior 
of crystalline 4 He at temperatures T < 200 mK. For 
the lack of any definitive experiment, we assume that 
the glass in 4 He is formed due to dislocations. In this 
Letter we would like to decouple the discussion of the 



observed features in specific heat and torsional oscilla- 
tor from the supersolidity. We will focus on two critical 
features reported to date by experiments: (1) the small 
entropy that is seen experimentally near the transition 
temperature. The observed entropy is at almost two or- 
ders of magnitude smaller than the expected entropy of 
a 1% supersolid fraction. (2) The observed linear specific 
heat in the bosonic crystal of 4 Hc, consistent with a glass. 

Entropy analysis. We focus on the specific heat mea- 
surements on 4 He by Clark and Chan [l4| which indi- 
cate a departure from the conventional T 3 specific heat 
behavior expected at low temperatures. Given the data, 
we searched for evidence of a thermodynamic phase tran- 
sition to a supersolid phase, assuming a 1% concentration 
of condensate. The observed features in the specific heat 
occur in the same temperature range where a change in 
the period of the torsional oscillator led to the specu- 
lation of a transition to a supersolid state (see Fig. HJ. 
A low-temperature linear term in the specific heat can 
be discerned. We note that a linear term in the specific 
heat of 3 He and 4 He crystals was observed more than 40 
years ago by Heltemes and Swenson and Franck , 
but not in later measurements by Edwards and Pandorf 
[l7j . A theoretical explanation of the earlier results was 
given in terms of thermal vibrations of pinned dislocation 
segments based on the theory of Granato [Isj . 

We now perform a general thermodynamic analysis 
of the measured specific heat data by Clark and Chan 
(see Figs. Q and |3J). This is especially attractive for 
the analysis as its result will not depend sensitively on 
the exact functional form of the specific heat C(T) at 
low temperatures, T < 100 mK. Most of the change 



C(T') 



which we find from the 



in the entropy, J dT' 
data, originates from the region of 100-400 mK. Accord- 
ing to the torsion oscillator measurements the supposed 
solid to supersolid transition occurs in the temperature 
range 160-320 mK. The low-T entropy differences of the 
760 ppm and 30 ppm samples relative to solid 4 He are 
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FIG. 1: S pec ific heat over temperature of solid 4 He by Clark 
and Chanfl4||. (Top): sample with 760 ppm 3 He at 38 bar. 
(Bottom): sample with 30 ppm 3 He at 39 bar. The dashed 
and solid lines are the respective low and high temperature 
fits, whose coefficients are presented in Table Q The insets 
show an enlarged temperature region. 



FIG. 2: The AC/T differences and the corresponding en- 
tropies AS = J* dT'AC/T' at low T. (Top): sample with 
760 ppm of 3 He at 38 bar for cooling and warming runs; the 
high-T (T > 200 mK) contribution of solid 4 He, C = BT 3 , 
with -B760ppm = 158 mJ/(K 4 mol) was subtracted. (Bottom): 
sample with 30 ppm of 3 He at 39 bar; the high-T contribution 
with -B30ppm = 155 mJ/(K 4 mol) was subtracted. 



roughly AS'(760ppm) ~ 0.06 - 0.4 mJ/(K mol) and 
AS , (30ppm) ~ 0.3 mJ/(K mol) respectively, between 80 
mK and 400 mK, sec Fig. El [3 ■ 

We assume, consistent with the experimental findings, 
that roughly 1% of the bulk sample transforms into a 
supcrfluid. Thus 99% of the sample is unaffected and 
behaves like a perfect Debye crystal below ~ 500 mK 
(T < 0£>/5O). Hence, the supersolid fraction would 
indicate an entropy of S = 46 mJ/(K mol). We now 
compare extracted entropy changes with what may be 
expected from a supersolid to solid transition (i) in a 
dilute gas, as suggested by a 1% supersolid component, 
and (ii) in a strongly correlated dense gas: (i) The Bose 
Einstein condensation (BEC) of an ideal gas gives val- 
ues of entropy changes which are an order of magnitude 
larger than those observed for the A transition of liquid 
4 Hc. Wc follow suite here and compare the ideal BEC 
with that of a possible supersolid transition. Assuming 
a noninteracting Bose-Einstein gas in three dimensions 
with a parabolic band, the specific heat of the BEC is 

°( t ) = tWM r{ ?/ Tc ^' 2 for T - Tc ' with c ^ the 

Riemann zeta function and gas constant R. The total 



entropy per mole in the condensed state taken between 
T = and T c is universally S B ec{T c ) = §f§7§-R ~ 5R. 
If we use a 1% molar concentration of supersolid fraction, 
we still get 5 B ec(T c ) « 0.05 R = 416 mJ/(K mol) [per 
mole of 4 He]. This is three orders of magnitude larger 
than the entropy changes seen in experiments (Fig. [2J • 
These conclusions do not change if a short-range repul- 
sive interatomic potential augments the kinetic energy in 
a dilute Bose gas. (ii) A lower bound can be found by 
scaling the measured S\ at the A-point of strongly corre- 
lated supcrfluid 4 He under pressure (Pa ~ 26 bar), which 
is S x = 4.6 J/(K mol) at T A = 1.8 K 0. As a gener- 
ous lower bound on the entropy at an actual supersolid 
transition, we may correct for the fact that the claimed 
transition temperature is 9 times smaller (~ 200 mK for 
30 ppm sample) 0] than T\. We scale down the expected 
1% entropy by a factor of 9. Then, we may expect to 
find an entropy of 5a ~ 5 mJ/(K mol). This estimate is 
about an order of magnitude bigger than the measured 
entropies. In ideal superfluids a nearly T c independent 
entropy S(T C ) results, even for different low-T scaling 
regimes, e.g., phonons, rotons. Wc emphasize that our 
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linear scaling of S(T C ) is a generous lower bound |23| . 

More sensitive measurements of the specific heat us- 
ing a new experimental setup by Chan et. al. (unpub- 
lished) are consistent with the deviations from a T 3 be- 
havior, however, the precise form at lowest temperatures, 
T < 80 mK is not settled yet. Clearly, the estimated en- 
tropy associated with an actual superfiuid transition is 
substantially larger than the measured AC/T integrated 
from ~ 80 mK to 500 mK for the case of A transition 
or weakly interacting BEC. The observed deficit of en- 
tropy is very hard to reconcile with the 1% fraction of 
superfiuid transition that is suggested to account for the 
torsional oscillator experiments. 

A null result of the signature of a A transition, the 
observation of hysteresis on cooling and warming and the 
dependence of the excess entropy on annealing points to 
a glass- like phenomenon of tunneling systems rather than 
that of a supersolid. 

Linear specific heat and sensitivity to 3 He impurities. 
We propose that the linear specific heat term 0,0,0 is 
due to Tunneling Systems (TS) in the 4 He crystal. Specif- 
ically, we assume that dislocation loops are small enough 
to be present in large numbers (~ 10 10 cm~ 2 ) to create 
the TS. We propose that fluctuations of nanoscale seg- 
ments of dislocation loops form the TS. A small addition 
of 3 He atoms to the 4 He crystal facilitates the creation of 
dislocations; hence, the linear specific heat term should 
increase with 3 He concentration, similar to the effect of 
small amounts of hydrogen in crystalline tantalum |2^j . 

The distribution of the characteristic energies of the 
TS is given by P(E) = PodE, where we follow the stan- 
dard discussion on the role of two-level systems in glasses 
pel Elf . It is commonly assumed that at low energy P 
is constant. We assume that Pq is only a function of the 
3 He concentration 713 and any contribution from intrinsic 
defects in 4 He is neglected. 3 Hc facilitates the creation 
of dislocations and hence Po should grow with the con- 
centration of 3 He. We further assume that per mole of 
4 He, Pa(ri3) = cNaTI^, with positive exponent v and co- 
efficient c, and Avogadro's number Na- A natural choice 
is v = 1, at least for low concentrations of 3 He. To keep 
the discussion general, we will not specify v. 

The specific heat of a single tunneling system is 
C TS (E,T) = k B (E/k B T) 2 [I jgg^g 7F , The aver- 
age over the distribution of the TS gives the total mo- 
lar specific heat, which at low temperatures is Cts = 
f™dE P(E) C T .s(E,T) ps ^-k 2 B P a T. The total spe- 
cific heat of 4 He is the sum of lattice, (C/ at ) and the TS 
(Cts) contributions. For a perfect Debye crystal, the 
molar C /at = ^R(T/Q D ) 3 at low T. In solid 4 He, the 
Debye temperature 9d ~ 28 K at P ps 40 bar. Thus at 
low temperatures (T < Od/50 ~ 0.5 K) the specific heat 
per mole of 4 He is C = AT + BT 3 , with A = ^k 2 B P and 
B = ^R/Q 3 D . For T<T* = #ppW(12^) ~ 
rig^ 2 , the linear term dominates over the lattice contribu- 



TABLE I: Summary of the linear and cubic coefficients of C = 
AT + BT 3 of solid 4 He with 3 He solute, as well as its Debye 
temperatures 6d. The A is from fits below ~ 120 mK, while 
the B and @ D are from fits between 200 mK < T < 500 mK. 
Rough estimates of the uncertainties are in parentheses. For 
30 ppm, A is close to zero as expected for this concentration. 

3 He P A B Qd 

ppm bar mJ/(K 2 mol) mJ/(K 4 mol) K 

30 39 0.09(4) 155(5) 23.2 

760 (cooling) 38 1.4(2) 158(5) 23.1 

tion. It follows that the TS model predicts that both the 
crossover temperature T* and the linear coefficient A in 
the specific heat will be very sensitive with respect to 3 He 
concentration. For example, for v = 1 it leads to a square 
root dependence on n 3 . This result also suggests that the 
effect of mass decoupling either vanishes or occurs at a 
much lower temperatures in samples with vanishing 3 He 
concentration. 

We revisited the data by Clark and Chan 0], see 
Fig. ^ an d assumed for our analysis linear and cubic 
terms in the specific heat. We find that indeed the data 
are consistent with a strong dependence on 3 He con- 
centration. For example, the linear term depends on 
3 He concentration, as shown in Table [I] The A coeffi- 
cient for 4 He reported by Heltemes and Swenson^] and 
Franck[Tr| varied between A ~ 2.5—8.8 mJ/(mol K 2 ) at a 
similar pressure of P ps 40 bar and with a Debye temper- 
ature ®£> ps 28 K.|23 Their A coefficients are somewhat 
bigger but of the same order as ours, which are in the 
range of A = 0.09 — 1.4 mJ/(K 2 mol), depending on 3 He 
concentration and history. 

The extracted low-temperature linear term coefficient 
A scales roughly linearly with the 3 He concentration. Al- 
lowing for different functional forms for the specific heat 
for various systems of interacting bosons does not lead 
to a significant change in our results. 

The "missing" moment of inertia and susceptibilities. 
The torsional oscillator experiments measure the suscep- 
tibility - they do not directly monitor the moment of 
inertia of the supersolid. As in any time translationally 
invariant system, the Fourier amplitude of the angular 
response of the torsion oscillator is 6(u>) — x(w)r ex4 (w), 
with x = Xi + *X2 an angular susceptibility and r ext 
the external torque. For the simple torsional oscillator, 
X _1 (o;,T) = [a - iuj"f osc - 1 osc lo 2 + g(uj,T)]. Here, I osc 
is the moment of inertia of the torsional oscillator, a is 
its restoring constant, 7 osc is the dissipative coefficient of 
the oscillator, and g(u),T) arises from the back action of 
solid 4 He. For an ideal normal solid with moment of iner- 
tia I ns , which rotates with the oscillator, the back action 
is g = —I ns (T)uj 2 . However, we do not need to impose 
this form. Current experiments measure the oscillator 
period 2it/ujo with ujq the real part of the solution of 
X~ 1 {to,T) = at fixed T. For example, a decrease in an 
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effective dissipative component (—i^fgi aS s^>) m g(u,T) as 
the temperature is lowered (wherein a liquid component 
transforms into a solid glass) will also lead to a shorter 
rotation period. The decrease in the rotation period only 
implies a crossover in x (and a constraint on g). As the 
real and imaginary parts of x are related by the Kramers- 
Kronig relations, an enhanced decrease in yi {y>, T) often 
appears with a pronounced peak in X2 ■ A nonvan- 
ishing xi at finite frequency mandates dissipation. 

The microscopic deformation of the glass to the applied 
torsion might, similar to suggestions |25[ concerning tor- 
sional oscillator results on granular media [2(|, take the 
form of elastic, plastic, and fracture processes (including 
those of internal avalanches). A gradual change can pro- 
ceed through dislocation glide in the slip plane [13] ■ All 
these possibilities need to be addressed experimentally 
before we can have a definitive microscopic picture of the 
possible glass state. 

One way to differentiate between glassy effects and a 
true thermodynamic p s is to measure the frequency de- 
pendence of p s . For the dislocation TS that we propose, 
p s {T,uj) will be a nontrivial function of frequency which 
vanishes as a power of frequency ps(T — ► 0, u>) ~ u> a . 
On the other hand, for a true supersolid phase there is a 
finite limit p s (T -> 0,oj) ~ p s (0). 

Heat Pulse Experiment. A heat pulse experiment 
may test for the existence and relevance of the TS for the 
thermal properties [2^|. In the absence of TS, the heat 
pulse would trigger quick equilibration of the energy with 
the phonon bath. Hence, on a very short time scale the 
temperature of a sample will reach a steady state value. 
By contrast, for a crystal with TS, the heat deposited in 
the crystal will first be absorbed by the phonon bath and 
later will cause the re-population of the TS as a result of 
the energy transfer from phonon bath to TS. 

Conclusions. We find that the measured entropy ex- 
cess is several orders of magnitude smaller than the en- 
tropy expected from the BEC or A transition of 1% su- 
perfluid fraction. While the absence of the entropy re- 
leased at claimed supersolid transition is puzzling, it is 
consistent with a dislocation induced glassy state in 4 He 
crystals. This hypothesis predicts that the linear term 
in the specific heat increases with 3 He concentration. 
Heat pulse, heat transport and ultrasound measurements 
would be helpful in identifying tunneling systems if they 
indeed exist in crystalline 4 He. The increase in the os- 
cillator frequency at low temperatures is consistent with 
this interpretation. 
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